home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 1998 November / IRIX 6.5.2 Base Documentation November 1998.img / usr / share / catman / p_man / cat3 / complib / zsytri.z / zsytri
Text File  |  1998-10-30  |  3KB  |  133 lines

  1.  
  2.  
  3.  
  4. ZZZZSSSSYYYYTTTTRRRRIIII((((3333FFFF))))                                                          ZZZZSSSSYYYYTTTTRRRRIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZSYTRI - compute the inverse of a complex symmetric indefinite matrix A
  10.      using the factorization A = U*D*U**T or A = L*D*L**T computed by ZSYTRF
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE ZSYTRI( UPLO, N, A, LDA, IPIV, WORK, INFO )
  14.  
  15.          CHARACTER      UPLO
  16.  
  17.          INTEGER        INFO, LDA, N
  18.  
  19.          INTEGER        IPIV( * )
  20.  
  21.          COMPLEX*16     A( LDA, * ), WORK( * )
  22.  
  23. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  24.      ZSYTRI computes the inverse of a complex symmetric indefinite matrix A
  25.      using the factorization A = U*D*U**T or A = L*D*L**T computed by ZSYTRF.
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      UPLO    (input) CHARACTER*1
  30.              Specifies whether the details of the factorization are stored as
  31.              an upper or lower triangular matrix.  = 'U':  Upper triangular,
  32.              form is A = U*D*U**T;
  33.              = 'L':  Lower triangular, form is A = L*D*L**T.
  34.  
  35.      N       (input) INTEGER
  36.              The order of the matrix A.  N >= 0.
  37.  
  38.      A       (input/output) COMPLEX*16 array, dimension (LDA,N)
  39.              On entry, the block diagonal matrix D and the multipliers used to
  40.              obtain the factor U or L as computed by ZSYTRF.
  41.  
  42.              On exit, if INFO = 0, the (symmetric) inverse of the original
  43.              matrix.  If UPLO = 'U', the upper triangular part of the inverse
  44.              is formed and the part of A below the diagonal is not referenced;
  45.              if UPLO = 'L' the lower triangular part of the inverse is formed
  46.              and the part of A above the diagonal is not referenced.
  47.  
  48.      LDA     (input) INTEGER
  49.              The leading dimension of the array A.  LDA >= max(1,N).
  50.  
  51.      IPIV    (input) INTEGER array, dimension (N)
  52.              Details of the interchanges and the block structure of D as
  53.              determined by ZSYTRF.
  54.  
  55.      WORK    (workspace) COMPLEX*16 array, dimension (2*N)
  56.  
  57.      INFO    (output) INTEGER
  58.              = 0: successful exit
  59.              < 0: if INFO = -i, the i-th argument had an illegal value
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ZZZZSSSSYYYYTTTTRRRRIIII((((3333FFFF))))                                                          ZZZZSSSSYYYYTTTTRRRRIIII((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              > 0: if INFO = i, D(i,i) = 0; the matrix is singular and its
  75.              inverse could not be computed.
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.